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ABSTRACT 

We report the results of a direct imaging survey of A- and F-type main sequence stars search¬ 
ing for giant planets. A/F stars are often the targets of surveys, as they are thought to have 
more massive giant planets relative to solar-type stars. However, most imaging is only sensi¬ 
tive to orbital separations > 30 AU, where it has been demonstrated that giant planets are rare. 
In this survey, we take advantage of the high-contrast capabilities of the Apodizing Phase 
Plate coronagraph on NACO at the Very Large Telescope. Combined with optimized princi¬ 
pal component analysis post-processing, we are sensitive to planetary-mass companions (2 
to 12 Mjup) at Solar System scales (<30 AU). We obtained data on 13 stars in L'-band and 
detected one new companion as part of this survey; an M6.0 + 0.5 dwarf companion around 
HD 984. We re-detect low-mass companions around HD 12894 and HD 20385, both reported 
shortly after the completion of this survey. We use Monte Carlo simulations to determine 
new constraints on the low-mass (<80 Mjup) companion frequency, as a function of mass and 
separation. Assuming solar-type planet mass and separation distributions, normalized to the 
planet frequency appropriate for A-stars, and the observed companion mass-ratio distribution 
for stellar companions extrapolated to planetary masses, we derive a truncation radius for the 
planetary mass companion surface density of <135 AU at 95% confidence. 
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1 INTRODUCTION 

Stellar properties are an important metric in the search for plan¬ 
ets, as they guide the target selection for detection surveys. In 
particular, stellar mass and metallicity are significant quantities in 
determining both the formation and evolution of stars and plan¬ 
ets (Johnson et al. 2010). Several radial velocity (RV) studies have 
shown that the giant planet frequency increases with stellar metal¬ 
licity (Santos et al. 2004; Fischer & Valenti 2005). The giant planet 
population as a function of stellar mass, however, is not consis¬ 
tent between different planet detection techniques (Quanz et al. 
2012; Clanton & Gaudi 2014; Vigan et al. 2012). While progress 
has been made in linking the RV and microlensing populations 
(Clanton & Gaudi 2014), this is a challenging problem involving 
the synthesis of different biases and parameter spaces covered by 
all the detection techniques. 

Gas giant planets (> 1 Mj^p) are the only directly imaged plan¬ 
ets thus far, due to their increased self-luminous thermal emission 

* Based on observations collected at the European Organization for Astro¬ 
nomical Research in the Southern Hemisphere, Chile, ESO under program 
numbers 089.C-0617(A), 089.C-0149(A). 


and decreased contrast at infrared wavelengths with the star. Planet 
populations derived from RV surveys are often extrapolated to 
larger orbital separations to analyze the frequency of giant planets 
in direct imaging surveys (e.g. Lafreniere et al. 2007; Biller et al. 
2013). 

Planet formation scenarios (Alibert et al. 2011) and simula¬ 
tions extrapolating RV planet populations (Crepp & Johnson 2011; 
Johnson et al. 2007) suggest that massive stars (>1.3 Mq) are the 
most favorable targets for directly imaging planets, since they 
have proportionally more material to form giant planets. Indeed 
many directly imaged planetary mass companions have been found 
around A or F stars; HR8799 bcde (Marois et al. 2008, 2010), f3 
Pic b (Lagrange et al. 2009, 2010), HD 95086 b (Rameau et al. 
2013a,c), HD 106906 b (Bailey et al. 2014). The detection of the 
HR8799 planets was the result of the Vigan et al. (2012) Interna¬ 
tional Deep Planet Survey. Most surveys, however, have yielded 
null results (Desidera et al. 2015; Chauvin et al. 2015; Janson et al. 
2013; Rameau et al. 2013b; Biller et al. 2013; Chauvin et al. 2010; 
Heinze et al. 2010; Lafreniere et al. 2007; Kasper et al. 2007). 
These null results are likely due to the lack of contrast at small or¬ 
bital separations, where most planets are expected to be found. Typ¬ 
ical detection limits for these surveys are 5-20 Mjup for > 30 AU. 
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Planets are rare at large orbital separations (Chauvin et al. 2010; 
Lafreniere et al. 2007; Nielsen & Close 2010) but at Solar System 
scales (< 30 AU), stars are largely unexplored. 

The main limitations for direct imaging are stellar “speckles” 
which can appear brighter than a companion (Hinkley et al. 2009). 
Coronagraphs are used in order to reach smaller angular separations 
around stars. They reduce the diffraction due to scattered stellar 
light in the telescope optics but at a cost of reduced throughput. The 
Apodizing Phase Plate (APP; Kenworthy et al. 2010; Quanz et al. 
2010, 2013) coronagraph suppresses the diffraction in a 180° wedge 
around a star, increasing the chances of detecting a very close-in 
companion. Several studies have demonstrated the APP’s capabil¬ 
ity of reaching < 30 AU (Meshkat et al. 2015; Ken worthy et al. 
2013; Quanz etal. 2011). 

We aim to probe down to Solar System scales (< 30 AU) 
around 13 A- and F-type main sequence stars in order to detect 
giant planets as well as set constraints on the planet frequency. 
We use the APP coronagraph on NACO at the Very Large Tele¬ 
scope (VLT), the L'-band filter, and optimized Principal Compo¬ 
nent Analysis (PCA) to achieve deep sensitivity limits (2 to 10 Mj^p 
at < 30 AU). 

In Section 2 we describe our target selection process, the coro- 
nagraphic observations, our data reduction method and how we de¬ 
termine the sensitivity of our data. In Section 3 we discuss the sen¬ 
sitivity achieved, our new detection of an M6.0±0.5 dwarf compan¬ 
ion to HD 984, and our re-detection of companions to HD 12894 
and HD 20385. We run Monte Carlo simulations to determine the 
probability distribution of our results, in order to compare different 
planet population models for A-type stars. Our conclusions are in 
Section 4. 


2 OBSERVATIONS AND DATA REDUCTION 

Our sample was carefully selected to derive the best possible con¬ 
straints on the frequency of giant exoplanets on Solar System 
scales: nearby, young, and massive main sequence stars. Young 
planets are still warm from their contraction (Spiegel & Burrows 
2012). By converting gravitational energy into luminosity, they are 
bright in the infrared. However, determining the age of a main- 
sequence star can be extremely challenging. One way to deal with 
this difficulty is to only select targets which are members of nearby 
associations with well established ages. If they are all bona-fide 
members of the group, we can assume the stars are of a similar 
age. Except for one', our targets are all members of nearby young 
moving groups or associations: p Pic Moving Group (23 ± 3 Myr; 
Mamajek & Bell 2014), Tuc-Hor Association (40 Myr; Kraus et al. 
2014), AB Dor Association (125 ± 15 Myr; Barenfeld et al. 2013). 
Nearby stars allow us to search for companions at smaller physi¬ 
cal separations. We aim to reach planet sensitivity on Solar Sys¬ 
tem scales, where we expect more planets to reside (Chauvin et al. 
2010; Lafreniere et al. 2007; Nielsen & Close 2010). Thus, we 
have selected only stars which are less than 66 pc away. 

At the time of selection, most^ of the targets were known to 
be single stars (Mason et al. 2011; Pourbaix et al. 2009) and not in 
the denser nucleus of their association. However, shortly after our 

* At the time of our observations, HD 984, was believed to be a 30 Myr 
member of Columba association. Based on our detection of a low-mass 
stellar companion to HD 984 and independent isochrone fitting, we estimate 
the age of HD 984 to be 115 ± 85 Myr (Meshkat et al. accepted). 

^ HD 20385 was known to have a wide binary 12" away. 


Table 1. Overview of stellar values used for each target. 


Target 

Mass (Mo) 

L' mag 

Spectral type 

Distance (pc) 

Age (Myr) 

HD 203 

1.40 

5.2 

F3V 

39.4 ±0.6 

23 

HD 12894 

1.39 

5.5 

F4V 

47.8 ± 1.0 

40 

HD 25457 

1.21 

4.3 

F6V 

18.8±0.1 

125 

HD 35114 

1.16 

6.2 

F6V 

8.3 ±0.9 

40 

HD 20385 

1.13 

6.4 

F6V 

49.2 ± 1.5 

40 

HD 102647 

1.9 

1.9 

A3Va 

11.0±0.1 

40 

HD 984 

1.18 

6.0 

F7V 

47.1 ± 1.4 

115 ±85 

HD 13246 

1.18 

6.2 

F7V 

44.2±0.9 

40 

HD 40216 

1.24 

6.2 

F7V 

54.4 ± 1.3 

40 

HD 30051 

1.38 

6.0 

F2/3IV/V 

63.6±4.2 

40 

HD 25953 

1.16 

6.6 

F5 

55.2 ±2.9 

125 

HD 96819 

2.09 

5.2 

AIV 

55.6 ± 1.7 

23 

HD 123058 

1.30 

6.7 

F4V 

64.1±3.5 

40 


Distances are extracted from pai'allaxes in the Hipparcos catalog (van Leeuwen 2007). The L'-band 
mag is converted from A!^-band mag in the 2MASS survey (Cutri et al. 2003) to using Cox (2000). 

The masses are from Casagrande et al. (201 1) and Chen et al. (2014). All ages are based on 
membership in nearby young moving groups or associations, taken from Mamajek & Bell (2014); 
Ki'aus et al. (2014); Luhman et al. (2005); Barenfeld et al. (2013), except for HD 984 which we 
compute in Meshkat et al. accepted. The bottom six targets were only observed in one APP 
hemisphere. 

survey was completed, a companion was discovered around one 
of our targets, HD 12894, by Biller et al. (2013) and Rameau et al. 
(2013b). Another target, HD 20385, was found to have a compan¬ 
ion shortly after our data were acquired (Hartkopf et al. 2012). 

2.1 Observations at the VLT 

Data were obtained for 13 targets from 2011 to 2013 (088.C- 
0806(B), 089.C-0617(A) PI: Sascha Quanz) at the Very Large Tele¬ 
scope (VLT)/UT4 with NACO (Lenzen et al. 2003; Rousset et al. 
2003) and the APP coronagraph (Kenworthy et al. 2010). The vari¬ 
ation in the observing time for each target depends on the observ¬ 
ing conditions on the night; if the observing conditions fell below 
a threshold during the night, the data acquisition was cancelled. 
Data were obtained with the L27 camera, in the L'-band filter (A 
= 3.80yum and A/i = 0.62^m) and the NB 4.05 filter (A - 4.051;rm 
and AT = 0.02^m) depending on the star’s L'-band magnitude. The 
visible wavefront sensor was used with each target star as its own 
natural guide star. We observed in pupil tracking mode to perform 
Angular Differential Imaging (ADI; Marois et al. 2006). We inten¬ 
tionally saturated the point spread function (PSF) core (on average 
out to ~0'.'08) to increase the signal-to-noise (S/N) from potential 
companions in each exposure. Unsaturated data were also obtained 
to calibrate photometry relative to the central star. 

The APP generates a dark D-shaped wedge on one half of a 
target. Excess scattered light is increased on the other side of the 
target, which is not used in the data analysis. Two datasets were 
obtained with different initial position angles (PA.) for full 360° 
coverage around the target star. Data were obtained in cube mode. 
Table 1 lists the stellar properties for each of our targets. Table 2 
lists the observing conditions for all the data obtained. 13 targets 
were observed in at least one APP hemisphere. 

2.2 Data Reduction 

A dither pattern on the detector was used to subtract sky back¬ 
ground and detector systematics from the raw data, as detailed in 
Kenworthy et al. (2013). Subtracted data cubes are centroided and 
averaged over. The two APP hemispheres obtained for each target 
must be processed separately, since they were observed on differ¬ 
ent nights and thus have different speckle noise patterns. Optimized 
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Table 2. Observing Log for NACOA^LT 088.C-0806(B) and 089.C-0617(A) 


Target 

Observation dates UT 
(Hem 1, Hem 2) 

Number of data cubes 

Total integration time (min) 

On-sky rotation (°) 

Average DIMM seeing (") 

HD 203 

2011 Oct 12, 2011 Nov 07 

154, 64 

76.8, 58.0 

46.65, 60.43 

0.90, 1.19 

HD 12894 

2011 Dec 10, 2011 Dec 24 

65,55 

55.6, 47.1 

31.97, 27.92 

0.78, 1.23 

HD 25457 

2011 Dec 11, 2011 Dec 21 

71,65 

43.4, 39.8 

28.63, 25.92 

1.44, 0.65 

HD 35114 

2011 Dec 13, 2012 Jan 02 

35,35 

30.3, 30.3 

31.98, 28.20 

0.75, 1.02 

HD 20385 

2011 Dec 21, 2012 Jan 08 

47, 35 

40.4, 30.3 

30.54, 24.61 

1.82, 1.25 

HD 102647 

2012 Jun 01, 2013 Apr 26 

47,56 

47.3, 56.2 

20.34, 24.71 

1.61,0.76 

HD 984 

2012 Jul 18, 2012 Jul 20 

59, 65 

54.5, 59.9 

47.41, 42.46 

0.65, 0.86 

HD 13246 

2011 Dec 07,- 

107,- 

90.9, - 

45.43, - 

0.99, - 

HD 40216 

2012 Jan 03, - 

11,- 

10.1,- 

11.74,- 

0.88,- 

HD 30051 

2012 Jan 07, - 

56,- 

48.0, - 

2.7,- 

1.68,- 

HD 25953 

2012 Jan 13,- 

79,- 

67.3,- 

42.76, - 

0.77, - 

HD 96819 

2012 May 02,- 

49,- 

47.3, - 

105.46,- 

0.64, - 

HD 123058 

2012 May 21,- 

62,- 

63.5,- 

29.44, - 

0.77, - 


Data are in chronological order based on first Hemisphere observed. The last six targets were only observed in one APP hemisphere. Targets are listed in the 

same order fQmo 1 


PCA was run on both of the APP hemispheres independently for 
each target, following Meshkat et al. (2014). This involves creating 
a linear combination of principal components (PCs) from the data 
itself in order to model and subtract away the stellar diffraction. 
Only the 180° D-shaped dark hemisphere was used in the PCA 
analysis. We fixed the number of PCs at approximately 10% the 
number of input frames, as this yields the optimal PSF subtraction 
close (< T.'O) to the star. We searched for point sources using this 
method for all 13 of our targets, despite not having full 360° cover¬ 
age for 6 of them. 

We injected fake planets into our data (before PCA process¬ 
ing) in order to determine the 5cr sensitivity limit for each target. 
Unsaturated data of the star was used to inject the fake planets. 
We scaled the unsaturated data to the same exposure as the sat¬ 
urated data. The star was added to the data with a contrast of 5 
to 12 mag in steps of 1 mag and from O'.'18 to T.'36 in steps of 
0.13. The outer radius limit was chosen because the field-of-view 
(FOV) of the APP is limited to only the upper quarter of the detec¬ 
tor (Kenworthy et al. 2010). The planet injected data was processed 
with PCA, de-rotated, and averaged over for the final image with 
North facing up. 

The final image was smoothed by a T/D aperture, in order to 
remove features which are not the expected planet size or shape 
(Amara & Quanz 2012; Bailey et al. 2013). We define the S/N of 
the injected planet to be the value of a single pixel at the location of 
the planet divided by the root mean square (rms) of a ring around 
the star at the angular separation of the planet, excluding the planet 
itself. Only the statistically independent pixels (one smoothing ker¬ 
nel apart) were used to compute the rms. Figure 1 shows the 5cr 
contrast curves for all 7 targets with two APP hemispheres and one 
target (HD 96819) with nearly full sky coverage in one APP hemi¬ 
sphere^. For targets with two APP hemispheres, fake planets were 
added at a fixed P.A. in each hemisphere. The average S/N of the 
injected fake planets in each hemisphere is used (at the same sep¬ 
aration). In the overlapping region between the two hemispheres, 
the number of frames varies slightly. However, since these regions 


^ HD 96819 has on sky rotation of 105?46. Since we are in ADI mode, the 
180° APP “dark hole” region rotates on the sky and only a ~ 74?9 wedge is 
missing. 



Figure 1. 5cr contrast curves for the targets with full 360° APP coverage 
around the star and one target, HD 96819, with nearly 360° coverage. 


are small and the number of frames never varies by more than 20%, 
the impact on the contrast curves is small. 

On average, we achieved a contrast of 9 mag at 0'.'4 and 10 
to 11 mag at > (X'6. The decreased sensitivity at < Qt.'l around 
HD 20385 is an outlier compared to the other targets. One possible 
explanation is the extremely bright companion detected at ~ O'.'8 
(discussed in Section 3.2), which affects the PCA component de¬ 
termination. 

We used the COND evolutionary tracks (Baraffe et al. 2003) 
to convert the contrast curves to planet mass detection limits 
(Figure 2). We were sensitive to planetary mass objects (< 12 Mj^p) 
at different projected separations, depending on the target distance. 
The outer radius for the sensitivity curves were based on the lim¬ 
ited FOV of the APP (Kenworthy et al. 2010). Thus, while the sen¬ 
sitivity curves appear to flatten out, we cannot extend these curves 
beyond T.'5 since we were not sensitive completely around the star. 
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Figure 2. Detection limits in lupiter masses versus projected separation in 
AU for all targets with full APP coverage (and HD 96819). We are sensitive 
to planet mass companions (< llMj^p) for all of our targets. HD 102647 
is more sensitive at smaller projected separation because it is much closer 
than the rest of the targets (11 pc). We do not plot the sensitivity beyond 
I'.'S for each target, due to the limited FOV of the APP. 


3 RESULTS 

We detect one new M6.0 ± 0.5 dwarf companion to HD 984, and 
re-detect companions to HD 12894 and HD 20385. To estimate the 
astrometry, we first determined the centroid of the three compan¬ 
ions. In this way, we verified that our star was well centered in our 
data. Since the APP has an asymmetric PSF, this step is crucial. We 
then injected fake negative companions at the location of the com¬ 
panion to determine the photometry and astrometry with error bars. 
We also varied the flux of the fake negative companions to cancel 
out the companion flux, which in some cases varies up to 20% due 
to atmospheric fluctuations. 

We iteratively converged on the P.A., angular separation, and 
A magnitude by varying the position and contrast and taking a 
minimization over the XjD aperture at the location of the compan¬ 
ion. For the very bright companions to HD 12894 and HD 20385, 
we determined the photometry and astrometry from ADI alone 
rather than PCA. Minor variations in the brightness and position 
of the companion in each frame can lead to “striping” in the final 
PCA processed image. This striping occurs when PCA fits the re¬ 
maining flux around of the companion after the fake companion is 
subtracted, since we never perfectly subtract the companion in in¬ 
dividual frames due to seeing variability. Table 3 lists the properties 
of the companions we detect based on our APP data. The error on 
the PA. includes uncertainties from true North orientation, based 
on direct imaging observations 0.5°, Rameau et al. 2013b). 


3.1 HD 12894 

Rameau et al. (2013b) and Biller etal. (2013) reported the dis¬ 
covery of HD 12894 B. Rameau et al. (2013b) detected a AL' = 
2.7 ±0.1 mag point source 14 AU from the star. They concluded 
that it was likely bound to the star, based on a non-detection in 1999 
2MASS data, because a background source would be detected due 
to proper motion of HD 12894. Based on the contrast and an age 
of 30 Myr, they concluded that the companion was a 0.8Mq K6 
star. Biller et al. (2013) concluded it was a co-moving 0.46 ± 0.08 
Mg companion 15.7 ± 1.0 AU away from its star, based on contrast 
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Figure 3. Top: HD 12894 ADI processed image, with the companion at 
O'/SI. Middle: HD 20385 ADI processed image, with the companion at 
0’/87. The shape of the companion PSF is due to the APP PSF. Bottom: 
HD 984 PCA processed image with the companion O'/19 from star. 
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Table 3. Companion properties of our targets. 


Target 

Date 

Separation (arcsec) 

P.A. (deg) 

AL' (mag)) 

Mass (Mq) 

HD 12894 B 

2011 Nov 24 

0.31 ±0.01 

240.24 ± 1.27 

2.89 ±0.14 

0.53 ±0.04 

HD20385 B 

2011 Dec 21 

0.87 ±0.01 

118.67 ±0.49 

2.52 ±0.10 

0.33 ±0.03 

HD984 B 

2012 July 18 

0.19 ±0.02 

108.9 ±3.1 

6.0 ±0.2 

0.11 ±0.01 


of AH = 3.0 mag. The discovery was made with the VLT/NACO 
instrument in L'-band as well as the NICI instrument on Gemini 
South in //-band, respectively. 

At the time our data were acquired, the companion was not 
known. We detect the companion with a AL' = 2.89 ± 0.14 
mag at a PA of 240.24±1°27 and separation of Oi'S! ± Oi'dl. 
This corresponds to a projected separation of 14.8 ±0.8 AU, with 
d=47.8 ± 1.0 pc for the distance to the star (van Leeuwen 2007). 
Figure 3 shows the ADI processed image of the companion. Using 
the COND evolutionary models (Baraffe et al. 2003) and the age of 
the star (40 Myr: Kraus et al. 2014), we estimate the mass of the 
companion is 0.53 ±0.04 Mg. Our analysis of this companion is 
consistent with Biller et al. (2013). 

3.2 HD 20385 

At the time we proposed to observe this target, it was a known 
wide binary, with a companion TOK 78 B 12" away. Hartkopf et al. 
(2012) reported the discovery of a new, close companion around 
FID 20385 at 0'.'88, Al =3.5 mag. Ay - 5.2 mag. They estimated 
the companion has a period of 200 years. Due to its variable radial 
velocity, they suggested the companion could be two unresolved 
companions. 

We detect the companion in our data at 0'.'87 ± O'.'Ol with a 
P.A. of 118167 ± 0149 and contrast of AL' = 2.52 ± 0.10 mag. Us¬ 
ing the stellar distance of 49.2 ± 1.5 pc (van Leeuwen 2007), this 
companion is at a projected separation of 42.8 ±1.8 AU. Figure 3 
shows the ADI processed image of the companion to HD 20385. 
The companion’s PSF clearly shows the APP PSF structure, with 
the bright lobes smeared due to the rotation on the sky. The age 
of the system is 40 Myr based on membership in the Tuc-Hor As¬ 
sociation (Kraus et al. 2014). Using COND models (Baraffe et al. 
2003), we estimate the companion to have a mass of 0.33 ±0.03 
Me. 

3.3 HD 984 

As reported in Meshkat et al. {accepted), we detected a companion 
around HD 984 in our APP data, as well as archival direct imaging 
data. Figure 3 shows our PCA reduced image with 20 PCs. We con¬ 
firmed the companion is co-moving with HD 984 and determined 
that is an M6.0 ± 0.5 dwarf based on SINFONI integral field spec¬ 
troscopy (see Table 3). 

3.4 Monte Carlo Simulations 

We ran 10000 Monte Carlo simulations of the target stars in order to 
determine the probability distribution of detecting substellar com¬ 
panions in our data assuming power-law slopes for the mass and 
semi-major axis distributions, for both planets and brown dwarfs 
(BDs), following Reggiani et al. submitted. For the BD distribution, 
we adopted the stellar companion mass ratio distribution (CMRD) 


A star planet frequency 



Figure 4. Detection probability distributions for our survey, assuming an 
outer-radius cutoff of 140 AU for the planet separation distributions (Reg¬ 
giani et al., submitted). With our null result, we rule out the A-star frequency 
beyond a cutoff of 135 AU at the 95% confidence level. 


from Reggiani & Meyer (2013) and a log-normal separation distri¬ 
bution (De Rosa et al. 2014). We included all of our targets in these 
simulations, including those with only one APP hemisphere cov¬ 
erage. In order to account for the targets without full sky coverage 
in our simulations, we multiplied the overall planet frequency per 
target by the fraction of sky coverage achieved. 

For the planets, we adopted a model distribution for A-type 
primaries which assumes the mass and separation distributions 
measured by RV surveys for solar-type stars, and extrapolated at 
larger separations (Heinze et al. 2010, based on Camming et al. 
2008). According to Gumming et al. (2008), the planet mass and 
semi-major axis distributions are modeled with power laws of in¬ 
dex Q'=-1.31 and /3=-0.61 {dN ~ M-'-^^dM and dN ~ a-°-^^da). 
To account for the A-type primary planet frequency, we adopted 
the values corresponding to the median sensitivity achieved in the 
RV data presented in Bowler et al. (2010) and Johnson et al. (2010) 
(f=l 1 ± 2% for planets in the ranges 0.5-14 Mj„p and 0.1-3.0 AU). 

For both planets and BDs, we assume a random distribu¬ 
tion of inclinations and the eccentricity distribution given by 
Juric & Tremaine (2008). In each simulation, we assigned each tar¬ 
get a number of planets and BDs from a Poisson distribution, ac¬ 
cording to the average number of planets and BDs per star, cal¬ 
culated from the aforementioned distribution. We also varied the 
outer radius cutoff: 20, 30, 80, 100, 120, 130, and 140 AU. The in¬ 
troduction of an upper limit for the planet separation distribution 
has been suggested by the results of previous direct imaging sur¬ 
veys (Chauvin et al. 2010; Vigan et al. 2012). 

If a target turns out to have one or more companions in the 
simulation, we assigned each companion a mass and the orbital 
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Table 4. Probability of a null result in our A and F main sequence star 
survey. 


^cutoff 

P(0) [%] A-type star planet frequency 

20 AU 

80 

30 AU 

70 

80 AU 

25 

100 AU 

14 

120 AU 

8 

130 AU 

6 

140 AU 

4 


parameters (semi-major axis, eccentricity, inclination) randomly 
drawn from the assumed distributions. The mass is converted into 
apparent magnitude, given the distance and age of the star and as¬ 
suming the same family of evolutionary models as in Section 2.2 
(COND Baraffe et al. 2003). The semi-major axis was converted 
into a projected separation, given the eccentricity and inclination 
and taking into account the time spent on the orbit. If the combi¬ 
nation of brightness and separation lies above the contrast curve 
(Figure 1), then the companion is detectable. Thus, at the end of 
every simulation, we know how many companions are created and 
how many are detected. After 10000 simulations, we determined 
the average detection probability for our A and F main sequence 
star survey (Figure 4). Three companions were detected in this sur¬ 
vey, but none of them were sub-stellar (<80 Mj^p). Thus, the proba¬ 
bility of detecting 0 companions for each model is given in Table 4. 

Given these probabilities, our null result allows us to re¬ 
ject the A-type star model with a scaled up planet frequency for 
^cutof f > 135 AU, at 95% confidence. This null result is also con¬ 
sistent with previous surveys that found that high mass planets at 
large orbital separations are rare (Nielsen et al. 2013; Biller etal. 
2013; Desidera et al. 2015; Chauvin et al. 2015). We note that un¬ 
like larger surveys, our simulations include both the planetary-mass 
companions as well as an extrapolation of the BD companion mass 
ratio distribution. Since we expect a contribution from both popu¬ 
lations in the total number of detections, a null result is more con¬ 
straining than a planet population alone. We also adopt a higher 
planet frequency (as expected for A type stars) compared to the 
standard planet frequency measured by Gumming et al. (2008) for 
solar-type stars. If the expected number of detections is higher, a 
null detection result allows one to place more stringent constraints. 
This is why our survey of 13 targets places comparable constraints 
compared to larger surveys. 

Based on RV measurements of A stars, Bowler et al. (2010) 
suggest positive values for the power law indexes of the mass and 
separation distributions, with high confidence. According to our 
Monte Carlo simulations, if we assume positive power law indexes, 
the probability of a null result is less than 0.1%, regardless of the 
planet frequency or the radius cutoff assumed. As suggested by 
Vigan et al. (2012), the inconsistency of direct imaging survey re¬ 
sults with the distribution parameters from RV observations around 
A-stars (Bowler et al. 2010) suggests that different planet popula¬ 
tions are probed by RV measurements at small separations than 
direct imaging at wide separations. 


4 CONCLUSIONS 

We present the results from a survey of carefully selected A- and 
F-type main sequence stars, searching for exoplanet companions. 


We aim to put direct imaging constraints on the occurrence of sub- 
stellar companions as a function of stellar mass. We obtained data 
on thirteen nearby (d<65 pc), young (<125 Myr) targets with the 
APR coronagraph on NACOWLT. We are sensitive to planet masses 
(2 to 10 Mjup) on Solar System scales (<30 AU) for all but one of 
our targets. We detected a new M6.0 ± 0.5 dwarf companion to HD 
984 and confirm stellar companions to HD 12894 and HD 20385, 
discovered shortly after our survey data were acquired. Our pho¬ 
tometry and astrometry for these companions are consistent with 
the values reported at the time of their discovery (Biller et al. 2013; 
Hartkopf et al. 2012). We found zero false positives in our L'-band 
data, as all of our detected point sources were bona-fide compan¬ 
ions. We perform Monte Carlo simulations to determine the ex¬ 
pected probability of detecting low-mass companions in our survey, 
based on our sensitivity and assumed semi-major axis distributions. 
Our non-detection of substellar companions (< 80Mjup) allows us 
to rule out the A-star frequency model distribution for >135 AU, 
with 95% confidence. 
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